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1  Introduction 

The  effects  of  heat  transfer  on  engineering  ap¬ 
plications  is  a  fundamental  aspect  in  the  design  of 
power  and  propulsive  systems.  The  heat  exchange  be¬ 
tween  fluid  and  solid  parts  may  reduce  the  efficiency 
of  the  system  in  situations  where  temperature  peaks 
and  gradients  become  important.  This  is  the  case  of 
gas  turbines,  combustion  engines  or  cooling  of  tur¬ 
bine  blades.  However,  this  interaction  is  difficult  to 
model  numerically  not  only  due  to  the  different  time 
scales  governing  each  individual  problem,  but  also  due 
to  instability  issues  related  to  the  coupling  strategy 
(Bunker,  2007).  The  coupling  of  the  Navier-Stokes 
equations  with  the  heat  equation,  also  referred  as  con¬ 
jugate  heat  transfer,  usually  requires  the  existence  of 
two  different  solvers,  one  for  each  domain,  exchang¬ 
ing  information  at  specified  boundaries.  Partitioned  or 
staggered  algorithms  can  be  applied  to  control  the  ex¬ 
change  of  information  allowing  the  solvers  to  advance 
at  different  time-steps  with  the  most  suitable  formula¬ 
tion  for  each  physical  problem  (Farhat  and  Lesoinne, 
2000). 

An  interesting  application  to  evaluate  the  numeri¬ 
cal  algorithms  and  the  effects  of  the  thermal  coupling 
with  the  flow  dynamics  is  the  case  of  a  jet  flame  im¬ 
pinging  on  a  flat  plate.  In  this  configuration,  high  rates 
of  wall  heat  transfer  are  achieved  as  a  consequence  of 
the  enhancement  of  the  heat  transfer  coefficients  by  the 
continuous  impact  of  vortical  structures  on  the  plate 
(Mira  et  al.,  2013).  Impinging  jets  are  characterised 
by  the  nozzle-to-plate  distance  H/D  and  the  dynamics 
of  the  flow  is  highly  dependent  on  this  ratio  (Hadziabic 
and  Hanjalic,  2008).  Three  flow  regions  can  be  iden¬ 
tified  in  this  configuration:  the  free  jet  region  where 
flow  is  not  influenced  by  the  plate,  the  deflection  re¬ 
gion  and  the  wall  jet  region.  For  large  H/D  ratios,  the 
configuration  permits  the  development  of  a  free  jet  and 
the  impingement  leads  to  the  formation  of  a  heat  trans¬ 
fer  profile  with  a  maximum  located  at  the  stagnation 
point.  However,  for  small  H/D  ratios,  the  jet  core  may 
be  longer  than  the  H/D  distance  causing  the  jet  to  be 
laminar  at  the  plate  featuring  low  levels  of  tubulence 
and  a  minimum  heat  transfer  coefficients  at  the  stagna¬ 
tion  point.  As  the  flow  accelerates  parallel  to  the  plate, 


the  wall  shear  stresses  increase  leading  to  a  maximum 
in  the  heat  transfer  profile  at  the  wall  (Fuchs  and  Hal- 
lqvist,  2009). 

The  aim  of  the  present  work  is  to  investigate  the 
heat  transfer  effects  and  flow  dynamics  of  an  imping¬ 
ing  flame  with  low  nozzle-to-plate  distance  when  the 
solid  plate  is  considered  non-adiabatic  and  the  heat 
flux  through  it  is  solved  by  the  heat  conduction  equa¬ 
tion. 

2  Mathematical  formulation 

Governing  equations 

The  governing  equations  describing  the  reacting 
flow  field  correspond  to  the  low  Mach  number  approx¬ 
imation  of  the  Navier-Stokes  equations  along  with  the 
species  conservation  equations  in  the  context  of  large- 
eddy  simulation  (LES). 

The  flow  field  is  filtered  in  space  using  a  box  fil¬ 
ter  given  by  A  =  V1^3,  where  V  represents  the  cell 
volume  and  a  Favre  filter  operator  is  employed  to  ac¬ 
count  for  the  density  variations.  The  filtered  governing 
equations  for  multi-species  reacting  flows  read: 
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where  p  is  the  viscosity  of  the  mixture,  t(u)  = 
\{Yu  +  Ytu)  —  |V  •  ul  is  the  deviatoric  part  of 
the  rate  of  deformation  tensor  and  I  the  identity  ten¬ 
sor.  pth  is  the  thermodynamic  pressure,  which  for  the 
open  flow  case  is  set  to  constant  and  equal  to  the  ex¬ 
ternal  pressure,  a  is  the  thermal  expansion  coefficient, 
usually  a  ~  1/T  (Avila  et  al.,  2011),  k  is  the  thermal 
conductivity  of  the  mixture,  and  hk  and  cP:k  are  the 


enthalpy  and  the  specific  heat  of  each  species  k.  The 
production  rate  of  the  kth  species  Cjk  is  obtained  as  a 
summation  of  all  reactions  involving  species  k  and  as 
only  one  chemical  reaction  is  included,  the  production 
rate  equals  the  rate  of  progress  Cok  =  &).  The  index 
m  =  (1  to  AT)  denotes  the  individual  species,  while 
the  superscript  sgs  refers  to  the  subgrid  scale  terms 
coming  from  the  filtering  operation. 

The  variables  with  superscript  sgs  correspond  to 
the  unresolved  momentum  transport  rsgs  =  p(uu  — 
uu ),  the  unresolved  heat  flux  hsgs  =  p(uT  — 
uT)  and  the  unresolved  species  mass  flux  <frsk9s  = 
p(uYk  —  uYk)  respectively.  The  subgrid  scale  mo¬ 
mentum  transport  is  modelled  using  an  eddy  viscos¬ 
ity  vt  obtained  by  the  Wall- Adapting  Local  Eddy- 
viscosity  (WALE)  model  proposed  by  Nicoud  and 
Ducros  (1999).  The  heat  flux  in  the  subgrid  scale  hsgs 
as  well  as  the  subgrid  scale  species  mass  flux  &k9S 
are  both  modelled  using  a  gradient  diffusion  approach 
(Mira  et  al.,  2013).  The  solid  part  of  the  domain  will  be 
treated  by  the  conjugate  heat  transfer  (CHT)  approach 
for  which  the  heat  equation  is  solved. 

dT 

PcP-^=V-(kVT)+q  (5) 

Coupling  approach  for  the  LES-CHT  simulation 

The  coupling  of  the  conjugate  heat  transfer  of  the 
solid  domain  with  the  LES  solver  is  accomplished  in 
such  a  way  that  the  fluid  imposes  the  heat  flux  to  the 
CHT  solver,  while  the  solid  imposes  the  skin  temper¬ 
ature  Ts  to  the  fluid.  The  boundary  conditions  are 
imposed  over  the  interface  using  the  nature  of  the  in¬ 
formation  received  by  each  solver.  The  heat  equa¬ 
tion  requires  a  Neumann-type  boundary  condition  q 
=  qs,  while  the  Low-Mach  fluid  uses  a  Dirichlet-type 
boundary  condition  for  the  temperature  T=TS  (Jaiman 
et  al.,  2006).  The  coupling  of  the  solid  and  fluid 
meshes  is  shown  in  Fig.  1.  Solid  nodes  are  represented 
by  thick  dots  on  the  interface  of  the  fluid  domain  and 
no  mathing  is  required  for  the  interpolation  algorithm. 

An  adaptation  of  the  Parallel  Location  and  Ex¬ 
change  library  (PLE)  from  Code_Saturne  intended  to 
provide  mesh  or  particle-based  code  coupling  services 
is  employed  for  the  LES-CHT  coupling.  It  allows  an 
easy  communication  between  application  codes  writ¬ 
ten  in  C/C++,  Fortran  or  Python  that  are  running  and 
share  the  same  communicator.  This  communicator  is 
split  so  that  there  is  one  communicator  for  the  exe¬ 
cution  of  each  solver  (LES  and  CHT).  The  communi¬ 
cation  between  the  applications  takes  place  using  the 
intracommunicators  created  after  spliting  the  global 
communicator. 

Chemical  kinetics  and  combustion  model 

The  chemistry  governing  the  oxidation  of  methane 
is  taken  into  account  by  the  1-step  chemical  kinetic 
mechanism  with  equivalence  ratio  correction  from 
Mantel  et  al.  (1996).  The  Thickened  Flame  model 


Figure  1:  Interface  matching  between  fluid  and  solid 
meshes. 


from  Collin  et  al.  (2000)  is  employed  to  thicken  the 
flame  front  so  that  it  can  be  resolved  in  the  LES  mesh. 
A  dynamic  calculation  of  the  local  thickening  factor  F 
is  employed  to  correctly  compute  the  flame  front  in  the 
mesh.  The  dynamic  thickening  factor  approach  from 
Durand  et  al.  (2005)  is  followed  to  obtain  the  local 
thickening  factor  F : 

F  =  1  +  {Fmax  —  1)11  (6) 


where  D  is  a  sensor  detecting  the  presence  of  the  flame 
front  and  Fmax  is  the  maximum  thickening  factor  al¬ 
lowed  in  the  computation.  This  sensor  can  be  ex¬ 
pressed  as  a  function  of  a  reaction  progress  c  that  is 
defined  as  the  ratio  of  mass  fraction  over  stoichoimet- 
ric  mass  fraction  c  =  1  —  Yf/Yf,st: 

Q(c)  =  16(c(l-c))2  (7) 


Fmax  is  calculated  as: 
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where  N  is  the  total  number  of  nodes  to  reproduce  the 
gradient  (N=6  in  our  current  simulation)  and  Si  is  the 
flame  thickness. 


Numerical  methods 

The  LES  equations  are  solved  using  the  Variational 
Multiscale  Stablization  (VMS)  technique  (Houzeaux 
and  Principe,  2008)  by  means  of  the  multiphysics  code 
Alya.  The  discretization  of  the  Navier-Stokes  equa¬ 
tions  yields  a  coupled  algegraic  system  to  be  solved  at 
each  linearization  step  within  a  time  loop.  Algebraic 
solvers  to  solve  this  coupled  problem  are  not  robust 
enough;  the  system  is  therefore  split  to  solve  the  mo¬ 
mentum  and  continuity  equations  independently.  This 
is  achieved  by  applying  an  iterative  strategy,  namely 
the  Orthomin  (1)  method  for  the  Schur  complement 
of  the  pressure  (Houzeaux  et  al.,  2011).  At  each  lin¬ 
earization  step  it  is  necessary  to  solve  the  momentum 
equation  twice  and  the  continuity  equation  once.  The 
momentum  equation  is  solved  using  the  GMRES  or 
BICGSTAB  method  (diagonal  and  Gauss-Seidel  pre¬ 
conditioners  are  usually  efficient),  and  the  continu¬ 
ity  equation  is  solved  using  the  Deflated  Conjugate 


Gradient  method  (Lhner  et  al.,  2011)  together  with 
a  linelet  preconditioner  well- suited  for  boundary  ley- 
ers.  The  parallelization  of  the  solver  is  extensively  de¬ 
scribed  elsewhere  (Houzeaux  et  al.,  2009).  The  GM- 
RES  solver  is  also  employed  to  solve  for  the  enthalpy 
and  species  mass  fractions.  The  Gauss-Seidel  iterative 
method  is  employed  to  solve  the  species  mass  fraction 
until  the  targeted  convergence.  The  code  has  been  val¬ 
idated  in  separate  efforts  (Houzeaux  et  al.,  2009)  and 
no  validation  results  will  be  shown  here  for  the  sake  of 
brevity. 

3  Computational  cases 

Physical  problem 

The  computational  cases  considered  in  this  work 
correspond  to  a  methane  premixed  flame  with  stoi- 
choimetric  equivalance  ratio  impinging  on  a  flat  plate. 
The  nozzle-to-plate  distance  is  H/D  =  2  and  the  inlet 
velocity  is  set  as  30  m/s  using  a  top-hat  velocity  pro¬ 
file.  The  jet  inlet  diameter  is  D  =  0.01  m  and  the  length 
of  the  fluid  domain  parallel  to  the  plate  is  10H.  Part  of 
the  inlet  pipe  with  length  H  is  also  accounted  for  in  the 
simulation.  The  solid  part  is  represented  by  a  squared 
domain  of  length  10H  and  thickness  H.  A  sketch  rep¬ 
resenting  the  fluid  and  solid  domains  is  shown  in  Fig. 
2. 
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Figure  2:  Sketch  of  the  computational  domain  separating 
fluid  and  solid. 


The  solid  domain  is  treated  as  a  metal  body  with 
the  thermal  and  physical  properties  of  copper.  In  or¬ 
der  to  keep  the  synchronization  in  physical  time  be¬ 
tween  the  fluid  and  solid,  the  metal  conductivity  was 
increased  by  a  factor  of  100.  It  means  that  the  solid 
response  to  changes  in  the  fluid  domain  is  enhanced 
allowing  to  obtain  strong  variations  in  temperature 
across  the  solid  and  better  visualize  the  cooling  ef¬ 
fects  on  the  impinging  flame.  Real  properties  of  the 
solid  domain  are  currently  being  investigated.  The 
boundary  conditions  in  the  solid  domain  are  adiabatic 
conditions  on  the  boundaries  and  a  fixed  temperature 
(Dirichlet-type)  Tont=450K  on  the  outer  part  of  the 
solid  domain. 

Impinging  wall  treatment 

Three  different  heat  transfer  conditions  at  the  im¬ 
pinging  wall  are  analysed.  Two  simplifications  that 
avoid  the  need  of  modelling  the  solid  and  the  entire 


coupled  problem  are  compared  to  the  real  case  where 
fluid  and  solid  are  solved  simultaneously.  The  first 
case  assumes  an  adiabatic  wall,  while  the  second  case 
imposes  a  Robin  condition  for  the  temperature  along 
with  the  heat  flux  across  the  plate.  For  this  approach, 
the  Tw  represents  the  wall  temperature  set  as  450  K 
and  h  is  the  convection  coefficient  set  as  50  W/m2K. 
The  last  case  computes  the  heat  flux  using  the  conju¬ 
gate  heat  transfer  (CHT)  approach  and  imposes  a  skin 
temperature  on  the  fluid  side,  and  the  heat  flux  on  the 
solid  side  (Duchaine  et  al.,  2008). 

The  boundary  conditions  specified  at  the  imping¬ 
ing  wall  (see  Fig.  2)  for  each  case  are  summarised 
here: 


Adiabatic 


Robin  type 
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Conjugate  heat  transfer  ( CHT) 
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Numerical  resolution  assessment 

The  resolution  in  the  fluid  domain  was  assessed 
by  testing  three  computational  meshes  of  1.0,  3.5  and 
5.2  million  elements  respectively  in  cold  flow.  The 
results  among  the  three  cases  were  compared  for  the 
mean  values  resulting  in  similar  trends  for  the  last  two 
cases.  The  maximum  normalized  wall  distance  y+  at 
the  impinging  wall  was  25,  15  and  10  respectively.  As 
the  Reynolds  number  reduces  in  reacting  flows  due  to 
the  increase  of  viscosity  inducing  an  increase  of  the 
boundary  layer,  the  mesh  with  5.2  million  elements 
was  selected  to  run  all  the  simulations.  The  resolution 
of  the  reacting  shear  layer  was  examined  by  the  local 
thickening  factor  obtained  in  the  simulations.  A  maxi¬ 
mum  local  thickening  factor  of  12  was  observed  in  all 
the  simulations  with  the  finest  grid  meeting  the  condi¬ 
tions  required  by  the  Thickened  Flame  model  (Collin 
et  al.,  2000). 


4  Results 

Instantaneous  fields 

The  main  flow  characteristics  of  a  premixed  im¬ 
pinging  flame  with  different  type  of  boundary  condi¬ 
tions  at  the  solid  wall  can  be  seen  in  Fig.  3.  It  shows 
the  three  main  regions  of  an  impinging  flame  (Mira 
et  al.,  2013):  the  free  jet  region,  the  deflection  zone 
and  the  development  of  a  wall  jet  using  a  scalar  rep¬ 
resentation  in  the  fluid  domain.  While  the  adiabatic 
case  shows  a  distribution  of  temperature  at  the  level  of 
the  cold  fuel  and  combustion  products,  the  cases  with 
heat  losses  indicate  the  existence  of  a  themal  boundary 


layer  where  the  temperature  reduces  substantially  to¬ 
wards  the  wall.  The  temperature  for  the  non-adiabatic 
cases  remains  more  uniform  over  the  wall  and  almost 
a  constant  value  can  be  observed  when  the  conjugate 
heat  transfer  approach  is  employed.  This  is  related  to 
the  fast  heat  conduction  over  the  solid  that  cools  down 
the  wall.  The  heat  losses  influence  the  formation  and 
convection  of  the  vortex  rings  as  well  as  the  dynam¬ 
ics  of  the  flame  in  the  deflection  region.  The  size  of 
the  vortex  rings  and  the  convection  of  these  vortices 
in  the  radial  direction  seems  to  be  affected  by  the  heat 
losses.  In  addition,  the  deflection  zone  is  shortened 
and  widened  by  heat  losses  affecting  the  formation  of 
the  vortex  rings  that  are  less  attached  to  the  wall  in  the 
adiabatic  case.  The  flow  dynamics  in  Fig.  3  indicates 
the  formation  of  a  wall  jet  accelerating  faster  and  sepa¬ 
rating  later  from  the  wall  for  non-adiabatic  cases.  This 
will  be  further  examined  using  the  time-averaged  flow 
fields  in  the  next  subsection. 
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Figure  3:  Temperature  contour  plot  for  the  adiabatic  (top), 
Robin  condition  (middle)  and  CHT  cases  (bottom) 
at  the  same  time  instant  t  =  0.006  s. 

Temperature  profiles  along  the  impinging  wall  also 
provide  fundamental  insights  into  the  heat  transfer 
characteristics  of  the  flames  under  investigation  and 
are  shown  in  Fig.  4.  The  adiabatic  case  shows  some 
unphysical  overshoots  in  temperature  caused  by  im¬ 
posing  a  zero  flux  over  the  boundary  with  no-slip 
condition  for  the  velocity.  The  heat  released  by  the 
chemical  reactions  increases  the  internal  energy  of  the 
fluid  without  any  dissipation  mechanism.  Flame/wall 
interactions  take  place  in  this  region  and  the  use  of 
adiabatic  conditions  to  account  for  this  phenomenon 
should  be  avoided.  The  plots  also  show  the  temper¬ 
ature  levels  and  distribution  along  the  plate  where  an 
important  reduction  in  temperature  is  observed  for  the 
cases  where  heat  losses  are  included.  Note  that  despite 
the  region  around  the  jet  centreline  is  rather  similar 
for  all  the  cases,  the  temperature  profiles  away  from 


this  region  are  quite  different.  The  addition  of  heat 
losses  in  the  treatment  of  solid  walls  acts  flattening  the 
temperature  distribution  and  reducing  the  temperature 
variations  on  the  interface.  A  fast  cooling  of  the  solid 
domain  is  observed  for  the  case  where  the  conjugate 
heat  transfer  is  solved  that  is  partially  enhanced  by  the 
increased  conductivity. 
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Figure  4:  Temperature  profiles  along  the  impinging  wall  for 
the  adiabatic  (top),  Robin  condition  (middle)  and 
CHT  cases  (bottom)  respectively. 


The  mean  temperature  and  heat  flux  across  the  im¬ 
pinging  plate  can  be  obtained  by  integrating  these  vari¬ 
ables  over  the  plate  surface  as  follows: 
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(13) 


where  dS  represents  the  differential  of  area  in  the  out¬ 
ward  normal  direction  to  the  surface.  The  time  evolu¬ 
tion  of  the  plate  mean  temperature  and  heat  flux  pro- 


vide  infomation  of  the  cooling  characteristics  of  the 
different  configurations  and  are  presented  in  Figs.  5 
and  6  respectively. 


Figure  5:  Mean  temperature  of  the  impinging  plate  over 
time. 
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products  that  is  eventually  reduced  to  the  equilibrium 
condition  as  the  flame  develops.  In  particular,  a  heat 
flux  coming  from  the  solid  to  the  flame  exists  during 
the  initial  transient  of  the  simulation  that  is  rather  large 
when  the  CHT  is  considered. 

The  distribution  of  temperature  gradients  in  the 
normal  direction  to  the  wall  over  the  impinging  plate 
relates  the  effect  of  the  heat  transfer  between  fluid  and 
solid,  and  is  shown  in  Fig.  7.  A  probability  distribu¬ 
tion  of  the  gradients  magnitude  over  the  normalized 
surface  area  is  also  shown  as  a  bottom  plot  in  Fig.  7. 
It  associates  the  magnitude  of  the  gradients  to  the  por¬ 
tion  of  area  on  which  they  are  distributed  providing  in¬ 
formation  related  to  the  sensibility  of  the  solid  respect 
to  changes  in  flame  temperature.  For  adiabatic  condi¬ 
tions,  a  zero  temperature  gradient  should  be  observed. 
However,  regions  where  the  heat  release  is  large  at  the 
wall  and  temperature  overshoots  appear,  a  temperature 
gradient  exists  to  balance  the  heat  equation  to  satisfy 
the  adiabatic  condition.  The  bottom  plot  shows  that 
these  temperature  gradients  are  rather  large  and  occur 
at  particular  locations  of  the  wall.  On  the  other  hand, 
for  non-adiabatic  conditions,  the  magnitude  of  the  gra¬ 
dients  reduces  and  extends  over  a  larger  part  of  the  sur¬ 
face.  This  is  further  examined  with  the  time-averaged 
fields  in  the  next  subsection. 
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Figure  6:  Mean  heat  flux  across  the  impinging  plate  over 
time. 


The  transient  evolution  of  the  temperature  and  heat 
flux  over  the  plate  shows  the  response  of  the  solid 
plate  to  the  flame  heating.  The  wall  mean  tempera¬ 
ture  for  the  adiabatic  case  shows  a  negligible  variation 
over  time  caused  by  the  the  condition  of  zero  flux  so 
no  heat  losses  are  taken  into  account.  The  case  with 
Robin  boundary  condition  shows  a  decay  in  temper¬ 
ature  respect  to  the  initial  condition  and  converges  to 
a  value  of  around  1150  K.  This  is  the  fluid  tempera¬ 
ture  at  the  start  of  the  thermal  boundary  layer,  since 
the  solid  temperature  was  imposed  at  450K.  The  case 
with  the  conjugate  heat  transfer  exhibits  the  fastest  de¬ 
cay  in  mean  temperature  over  the  plate  with  a  mean 
value  around  500K.  The  response  of  the  solid  to  the 
flame  heating  is  achieved  by  the  heat  flux  across  the 
solid  and  this  is  shown  in  Fig.  6.  The  case  with  the 
CHT  is  subjected  to  the  largest  heat  flux  and  there¬ 
fore  to  the  lowest  skin  temperature,  which  is  partially 
caused  by  the  high  conductivity  imposed  on  the  solid 
material.  The  profiles  for  the  three  cases  indicate  a 
large  heat  flux  at  the  start  of  the  simulation  caused  by 
the  initialization  of  the  domain  with  hot  combustion 


Figure  7:  Instantaneous  temperature  gradient  (top)  and 
probability  density  function  of  temperature  gradi¬ 
ents  (bottom)  at  the  impinging  plate  at  time  t  = 
0.015  s. 


For  the  CHT  case,  the  distribution  of  temperature 
over  the  fluid  and  solid  domains  is  shown  in  Fig.  8 
for  a  given  time  instant  in  a  streamwise  cross  section. 
The  core  of  the  flame  around  the  centreline  creates  a 
central  region  of  cold  temperature  on  the  solid  side, 
while  a  hot  region  is  found  off  side  the  center  where 
the  secondary  vortices  are  formed  and  the  combustion 
heat  release  interacts  with  the  solid.  The  heat  is  dis¬ 
tributed  across  the  domain  and  the  temperature  distri¬ 
bution  along  the  centreline  is  shown  in  Fig.  9.  The 


temperature  rapidly  reduces  from  the  impinging  wall 
to  the  fixed  value  set  on  the  outer  part  of  the  solid  do¬ 
main  due  to  the  high  conductivity  imposed  on  the  solid 
material. 


Figure  8:  Temperature  contour  plot  for  the  fluid  and  solid 
domains  at  time  t  =  0.010  s. 
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Figure  9:  Temperature  across  the  center  of  the  fluid  and 
solid  domains  at  different  time  instants. 

The  heat  transfer  condition  also  affects  the  heat  re¬ 
lease  distribution  in  the  near-wall  region  as  seen  in 
Fig.  10.  In  the  free  jet  region,  the  heat  release  rate 
is  rather  similar  for  all  the  cases,  but  the  differences 
start  appearing  in  the  deflection  zone  where  the  vortex 
rings  are  formed.  The  adiabatic  condition  enhances 
the  chemical  reaction  rates  by  sustaining  a  high  tem¬ 
perature  region  that  promotes  the  release  of  heat.  Fig¬ 
ure  10  also  shows  that  the  adiabatic  case  undergoes 
an  enlargement  of  the  reacting  region  causing  the  dis¬ 
similar  temperature  distribution  observed  among  the 
cases.  The  case  with  the  conjugate  heat  transfer  devel¬ 
ops  a  thin  region  with  chemical  reactions  taking  place 
changing  substantially  the  heat  release  distribution. 

Time-averaged  fields 

The  time  integration  of  the  instantaneous  fields  is 
obtained  after  time- averaging  the  data  for  two  flow 
through  calculations  using  the  averaged  radial  veloc¬ 
ity.  The  time  integration  is  considered  after  the  fields 
are  statistically  independent  and  no  effect  of  the  ini¬ 
tial  conditions  exist.  The  time-averaged  fields  are  em¬ 
ployed  to  investigate  the  heat  transfer  effects  on  the 
formation  of  the  wall  jet  and  the  distribution  of  tem- 
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Figure  10:  Heat  release  contour  plot  for  the  adiabatic  (top), 
Robin  condition  (middle)  and  CHT  cases  (bot¬ 
tom). 

perature  gradients. 

An  important  phenomena  in  impinging  flows  is  the 
formation  of  a  wall  jet.  To  examine  the  effects  of  the 
wall  heat  transfer  on  the  formation  and  development  of 
wall  jets,  velocity  profiles  at  different  radial  locations 
are  shown  in  Fig.  11.  Despite  the  velocity  distribution 
is  rather  similar  among  the  cases  near  the  centreline, 
a  further  accelaration  of  the  flow  occurs  for  the  cases 
with  heat  losses,  specially  for  the  CHT  case.  The  flow 
with  the  conjugate  heat  transfer  has  the  largest  peak 
velocity  and  the  weakest  back  flow  satisfiying  the  con¬ 
servation  of  momentum.  Results  not  shown  here  indi¬ 
cate  that  the  time-averaged  shear  stresses  are  the  same 
for  all  the  cases  suggesting  that  the  differences  in  size 
between  the  vortex  rings  are  caused  by  the  the  dissim¬ 
ilar  distribution  of  temperature  gradients  on  the  wall. 
The  normalized  distribution  of  temperature  gradients 
on  the  impinging  wall  is  shown  in  Fig.  12.  The  nor¬ 
malization  is  done  with  the  same  values  as  Fig.  7,  so 
a  direct  comparison  in  both  magnitude  and  distribu¬ 
tion  can  be  made.  Figure  12  evidences  the  reduction 
in  magnitude  of  temperature  gradients  as  well  as  an  in¬ 
crease  in  the  area  over  which  the  temperature  gradients 
are  distributed  from  the  Robin  case  to  the  CHT  case. 

5  Conclusions  and  future  work 

The  study  presents  some  large-eddy  simulation  re¬ 
sults  of  a  premixed  impinging  jet  for  a  configuration  of 
nozzle-to-plate  distance  of  H/D  =  2.  The  results  pre¬ 
sented  are  focused  on  the  flow  behaviour  and  the  heat 
transfer  effects  near  the  impinging  plate  when  differ¬ 
ent  boundary  conditions  are  considered  and  the  conju¬ 
gate  heat  transfer  is  employed  on  the  solid  domain.  A 


Wall  jet  velocity  (m/s) 


Figure  1 1 :  Wall  jet  development  using  time- averaged  flow 
fields. 


Figure  12:  Normalized  distribution  of  mean  temperature 
gradients  at  the  impinging  wall. 


new  coupling  algorithm  for  parallel  communication  is 
used  to  exchange  data  between  the  fluid  and  the  solid 
domain  in  a  efficient  manner  and  without  reducing  the 
parallel  efficiency  of  the  individual  solvers.  The  re¬ 
sults  reveal  substantial  differences  when  the  boundary 
condition  at  the  impinging  wall  is  changed.  Adiabatic 
conditions  lead  to  a  solid  wall  at  high  temperature  with 
a  large  variation  in  the  radial  direction,  while  the  non- 
adiabatic  cases  show  a  more  uniform  temperature  dis¬ 
tribution  over  the  plate.  The  effect  of  heat  losses  not 
only  tends  to  reduce  the  magnitude  of  the  temperature 
gradients,  but  also  to  enlarge  the  area  where  the  tem¬ 
perature  gradients  exist.  This  affects  the  wall  jet  de¬ 
velopment  accelerating  the  flow  parallel  to  the  plate 
and  changing  the  size  of  the  vortex  rings,  and  the  dy¬ 
namics  of  the  flames  near  the  wall.  The  future  work 
is  focused  on  investigating  the  properties  of  the  solid 
domain  in  the  dynamics  of  the  flame  and  the  stability 
of  the  coupled  problem. 
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